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ABSTRACT 

We investigate the dynamics of magnetic fields in spiral galaxies by performing 3D 
Magnetohydrodynamics (MHD) simulations of galactic discs subject to a spiral poten- 
tial using cold gas, warm gas and a two phase mixture of both. Recent hydrodynamic 
simulations have demonstrated the formation of inter-arm spurs as well as spiral arm 
molecular clouds provided the ISM model includes a cold HI phase. We find that the 
main effect of adding a magnetic field to these calculations is to inhibit the forma- 
tion of structure in the disc. However, provided a cold phase is included, spurs and 
spiral arm clumps are still present if jS > 0.1 in the cold gas. A caveat to the two 
phase calculations though is that by assuming a uniform initial distribution, /3 > 10 
in the warm gas, emphasizing that models with more consistent initial conditions and 
thermodynamics are required. Our simulations with only warm gas do not show such 
structure, irrespective of the magnetic field strength. 

Furthermore, we find that the introduction of a cold HI phase naturally produces 
the observed degree of disorder in the magnetic field, which is again absent from 
simulations using only warm gas. Whilst the global magnetic field follows the large 
scale gas flow, the magnetic field also contains a substantial random component that is 
produced by the velocity dispersion induced in the cold gas during the passage through 
a spiral shock. Without any cold gas, the magnetic field in the warm phase remains 
relatively well ordered apart from becoming compressed in the spiral shocks. Our 
results provide a natural explanation for the observed high proportions of disordered 
magnetic field in spiral galaxies and we thus predict that the relative strengths of the 
random and ordered components of the magnetic field observed in spiral galaxies will 
depend on the dynamics of spiral shocks. 
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1 INTRODUCTION 

Observations of the magnetic field in galaxies indicate 
that the magnetic pressure is comparable in magni- 
tude to the thermal pr essure and interstellar turbulence 
(|Heiles &: Crutche3 l2005h. Consequently magnetic fields are 
expected to play an important part in the dynamics of the 
ISM, including spiral shocks and the formation of molecular 
clouds. 

The response of a gaseous disc to a spi- 
ral potential has been well examined both 
theoretically and numerically (RobertsI 



shocks, providin g the potential is of sufficient strength. 
[Roberts Yua^ (|l970l ) found that the strength of the 
shock decreases with increasing magnetic field strength, 
whilst the magnetic field strength is amplified in the spiral 
shock, the latter agreeing with observations at the time for 
galactic magnetic fields. 

There are numerous hydrodyna mic simulations of gas 
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into a quasi- stationary solution which contains spiral 



discuss the formation of substructure and location of the 
shock. Whilst 2D simulations of warm gas show the for- 
mation of small-scale spurs perpendicular to the spiral arms 
(jGittins k Clarkell2004l : IWada k Kodall2004ll . these features 
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appear to be suppressed in 3D simulations (Kim k Ostr ike3 
2006). On the other hand, Dobbs & Bonnell (2006.) do 
find more extensive substructu r e and spiral arm clouds 
(jDobbs. Bonnell k Pringi3 l2006l : iDobbs k Bonnelll l2007bh 
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in 3D simulations, but only in models where the gas tem- 
perature is ^ 1000 K. There are comparatively fewer MHD 
calculations, largely due to the relatively recent addition of 
magnetic fields into n umerical codes, and the li mitations of 
numerical resolution. IShettv Ostriker (2006) investigate 
the formation of spurs using 2D MHD grid-based calcula- 
tions of the warm ISM with self gravit y, concluding that self 
gravi ty is necessary for spurs to evolve. [Comez Coxl (|200l 
l2004l l also perform 3D calculations with warm gas, but 
whilst both groups compare MHD with non-MHD results, 
neither investigate the effects of varying the field strength 
or ISM temperature. 

Recent observations now provide much more detailed 
measurements of the magnetic field strength and direction 
in t he diffuse (w arm) ISM (see reviews bv lBeck et al.lll996l 
and iBeckl l2007l l . Spiral magnetic arms appear to occur in 
all disc g alaxies, reg ardless of the presence of optical spi- 
ral arms (|Beckll2005l ). Typically the field in the inter-arm 
regions is several /xG, but may be 10 or more /xG in the 
spiral arms, and the total field contains random and or- 
dered (regular) components of comparable strength. The 
general consensus currently favours dynamo theory to ex- 
plain the orig in of these relatively strong magnetic fields 
(IParkedflQTlal lbl: iBalsara. Kim. Mac Low k M athews 200i; 
lBeckll2005l ). whereby the magnetic field is generated by tur- 
bulence in the IS M, and can produce or enhance magnetic 
spira l arm arms (IPanesar NelsonI [Tooi : iRohde Elstnerl 
1 199 8) . However in grand design galaxies where the gas expe- 
riences strong spiral shocks, such as M51, the ma gnetic spi- 
ral arms tend to be aligned with the dust lanes (jNeiningerl 
Il992l ) suggesting that the magnetic field is strongly related 
to the dynamics of the spiral shock. 

Magnetic fields have been long associated with molec- 
ular cl oud formatio n, originally through the Parker insta- 
bility (jParker 1 19661). and more recently via the magnet o- 



rotational ('Kim. Ostriker Stonell2003l : IPiontek Ostrikerl 
I2OO5I I instability. The Parker instability is expected to 
produce sinusoidal motions in the z direction, with the 
magnetic field channeling gas into dense concentrations 
in the plane of the disc. Molecular clouds formed in this 
way were originally believed to be magnetically supported, 
with lifetimes of order 10^ y ears (Zweibel & Kulsrud 1975; 
IShu. Adams Sz Lizanol [l987n . However the Parker instabil- 
ity is found to be relatively weak, and s ingularly insuffi- 
cient to induce molecular cloud formation (lElmegreenlll982l : 
'Kim. Rvu Jones 2001: "Kim. Ost riker Stonel '2002\ The 
magneto-rotational instability (MRI) has been p roposed as 
a fo rmation mechanism away from spiral arms (|Kim et al.l 
1 2 0031 ) . However the MRI generally takes a few orbits to ini- 
tialise, particularly if the field is t oroidal (e.g. lO's of or- 
bits in simulations by iNishikori et al.l 120061 ) , and is there- 
fore of less importance where spiral shocks occur. Mag- 
netic fields may be more relevant to GMC (giant molecular 
cloud) formation in conjunction with gravitational insta- 
biliti es (Elmegreen 1987; Kim k Ostriker 2001; Kim et al.1 
I2OO2I ) [a lso known as the magneto - Jeans-instability, MJI] or 
cooling (jKosihski k Hanasd [20071 ). 

In this paper we describe 3D MHD calculations of 
a galactic disc subject to a spiral perturbation. These 
are the first fully Lagrangian MHD calculations to model 
this problem, using the Smo othed Particle Magneto- 
hydrodynamics (SPMHD) code (jPrice k MonaghanI [2004bl . 



I2OO5I : IPrice k Batell2007l ). We compare the structure of the 
disc for a range of magnetic field strengths, assuming a single 
phase medium of cold or warm gas, or a two-phase medium. 
We also describe the strength and morphology of the mag- 
netic field, and relate these to observations. In this paper 
we focus on the effect of the spiral shock inducing structure 
in the gas rather than Parker or MRI instabilities, and we 
leave a discussion of results including self-gravity to a future 
paper. 



2 CALCULATIONS 

We use a 3D Smoothed Particle Magnetohydrodynamics 
code (SPMHD) for these c alculations, based on a binary tree 
code developed by Benz (|Benz et al.lll99Ql ). The code has 
been extensive ly modified, including the use of sink parti- 
cles (Bate 1995), and in particular, the addition o f magnetic 
fields (Price k Monaghan 2004b. al: [Price k Batel [2bo7). 

The SPMHD code solves the equations of magnetohy- 
drodynamics in the form 
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where p, , P, P and B are the fluid density, velocity, 
pressure, thermal energy per unit mass and magnetic flux 
respectively, = -^-\-v^ is the time derivative following 
the motion of a fluid particle and S^-^ is the stress tensor 
given by 



-PS'^ + — ( B'B 
po 



2 



(6) 



The numerical representation of these equations as 
su mmations over SP H parti cle s is discus sed in detail 
in iPrice k MonaghanI (|2004bl lal. l2005l ) and IPrice k Batl 
(|2007l ). In particular equations dH)-® use t he 'Euler po- 
tentials' represen tation for the magnetic field (|Price k Bate! 
l2007l : ISternlll97 0"). for which the Lagran gian evolution is zero 
(Equation [5]), corresponding to the advection of magnetic 
field lines in ideal MHD. The advantage of using this for- 
mulation is that the divergence constraint on the magnetic 
field is satisfied by construction. 

Dissipative terms in the form of an artificial viscosity 
and resistivity are added to the momentum and induction 
equations on order to correctly capture discontinuities in 
the flow (i.e. shocks and magnetic current sheets). Artifi- 
cial viscosity is used wi th the standard pa rameter s ay — 1 
and pv =2 ((MonaghanI 1 19971 : IPrice k MonaghanI I2OO5I ). as 
developed for the SPMHD code and artificial resistivity is 
added to th e Euler poten t ials' evolution with as — 1 as 
described in "Price & Batd (|2007l ). 

A galactic potential is added via an external driv- 
ing force fext, the exa ct form of which is described in 
iDobbs k Bonnelll (|2006l ). The main features of the spiral 
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component ffrom lCox G6mezll2Q02l ) are a 4 armed spiral 
pattern, of pattern speed 2 x 10^ rad yr~^ and pitch angle 
15^^. The amplitude of the spiral potential is 1.1 x 10^^ cm^ 
s~^. The total potential also includes components for the 
disc and halo. 

The code also includes individual time-steps, and vari- 
able smoothing lengths. The calculation of the smooth- 
ing lengths uses a recent improvement by which the den- 
sity and smoothing are so l ved i teratively for each particle 
(|Price Mo naghan 20Q4b'.'20Q7). The density and smooth- 
ing length are related according to 

where h is the smoothing length, p the density, m the mass 
of the particle and 77 is a dimensionless parameter set to 1.2 
in order that each particle has ^ 60 neighbours. This gives 
a mean value oi h ^ 30 pc initially (40 pc in the multi-phase 
simulations) using 4 million particles, such that h/H ^ 0.2. 
Although the scale height (H) of the cold gas decreases (see 
next section) with time, the density of the cold component 
increases considerably in the spiral arms, so h also decreases 
in the arms. Overall the z dimension tends to be marginally 
resolved in the spiral arms but not well resolved in the inter- 
arm regions. 



2.1 Initial distribution of gas 

We place particles in a torus with radii 5 kpc < r < 10 
kpc. The gas is initially distributed randomly and the ve- 
locities are assigned with circular orbits. Unlike previous 
results (Dobbs & Bonnell 2006), where the initial positions 
and velocities were extracted from a test calculation, the 
spiral perturbation emerges with time in these simulations. 
The velocities in the plane of the disc follow a rotation curve 
corresponding to the disc component of the potential 



-vo log[r'^ + Rl + {z/zqf 



(8) 



where Rc—l kpc, vq — 220 km s~^, and Zq—0.7 is a mea- 
sure of the disc scale height. This produces an essentially 
fiat rotation curve for the radii over which the particles are 
distributed. The velocities in the vertical direction are cho- 
sen from a Gaussian of mean 6 km s~^. The total mass 
in the torus is 10^ M©, giving an initial surface density 
of 4 Mopc"^, similar to that of HI near th e solar radius 
(|Wolfire. McKee. Hollenbach k TielensI l2003l ). In all these 
calculations we use 4 million particles. 

We perform single phase and two phase calculations. 
The gas is either cold (100 K) or warm (10^ K) in the sin- 
gle phase models. For the two phase calculations, we use 
equal proportions of warm and cold gas. Each phase con- 
tains 2 million particles placed in a uniform random distri- 
bution over the disc. As the simulation progresses the dis- 
tribution evolves in to cold clumps surr ounded by a warm 
tenuous medium (jPobbs Bonnellll2007bl V All calculations 
are isothermal, and we do not include self gravity or feedback 
from star formation. The assumption of an isothermal gas is 
not ideal, and for the two-phase simulations, we would ex- 
pect heating and cooling between the two components. Our 
analysis instead assumes that ther e is a reservoir of cold 
HI in the ISM as well as warm gas (jHeiles k Trolandlliooi : 



iGibson et aLll2006l ). and we consider the response of such a 
distribution to a spiral density wave. 

In the single phase simulations, we take an initial maxi- 
mum height of 150 kpc, whilst for the two phase simulations 
the initial maximum height is 400 pc. These values were cho- 
sen b ased on the observed scale heights for the Galaxy ([Co3 
I2OO5I ). for cold and warm gas. However assuming pressure 
equilibrium, and that the densit y distribution corresp ond- 



ing to Eqn. 8 falls off as z ( Bin nev Tremaini 
scale height expected from the logarithmic potential is: 



a c orresp o 
i^ ll98/1 ). 
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at a radius of 7.5 kpc, where Cs is the sound speed of the 
gas. This leads to a scale heights of approximately 33 and 
330 pc in the cold and warm phases respectively at a radius 
of 7.5 kpc. During the simulation, the scale height of the 
warm gas settles into an exponential distribution of scale 
height 300-340 pc. We note that without supernovae, a lower 
scale height for the wa rm gas is exp ected than implied by 
observations (- 130 pc ( McKeelllQQQl )). in which case a lower 
value of Zq may be appropriate (0.25-0.3). 

The scale height of the cold gas is initially maintained 
by the velocity dispersion applied in the z direction. How- 
ever the increase in density in the spiral arms and decay of 
the velocity dispersion means that the scale height of the 
cold gas is typically 20-100 pc at later times in both the sin- 
gle and two phase calculations. In the two phase simulations 
we did not assign a different initial height for the cold and 
warm components. However the results are not particularly 
sensitive to the initial scale height and do not change signif- 
icantly if we set the initial scale height of the cold gas to a 
lower value. 

The initial configuration of the magnetic field is 
toroidal, similar to previo us calculations (jShettv Ostriked 
|2006VG6mez & Cox"2004). A toroidal field can be described 
in terms of Euler potentials by 
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where — + y"^ + z^ and Q — cos~"^(z/r). The relative 
strength of the magnetic field is given by the mean plasma 
beta, defined as the ratio of gas to magnetic pressur4j 



/3 = 



BV2/io 



(12) 



where po is the average density of the disc. The Alfven 
speed is then va — B j ^fji^p^ — V^Cs/y/p. With 
prese nt observations of the diffuse (non-self-gravitating) 
ISM (|Heiles k Crutcherll2005l ) . the magnetic pressure is esti- 
mated to be 3 times larger than the thermal component. The 
CNM and WNM are both found to exhibit both magnetic 
field stren gths of ^ 6^ G and a th ermal pressure of ^ 10400 
cm"^ K ( Heiles k Trolandl I2OO5I 1 indicating that will be 
similar for each component. However in the two-phase sim- 
ulations described in this paper (Section 3.4), we assume 



The reader should note that in some papers (e.g. 
IShettv k Ostrikerl [ 20061 ). /3 is defined as half this ratio in 
which case va = Cg/^/p. 
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Model 


T(K) 




B (/iG) 


VA (km s ■*■) 


A 


100 


10^ 


3 X 10-^ 


1.2 X 10-3 


B 


100 


10 


0.09 


0.35 




iUU 


1 

i 


U.o 




D 


100 


0.1 


0.9 


3.5 


E 


10^ 


100 


0.3 


1.2 


F 


10^ 


1 


3 


12 



Table 1. Table showing the parameters used in the single phase 
calculations. The magnetic field strengths correspond to a uni- 
form initial ISM density of ~ 0.7 cm~^ . 

initial conditions of a uniform density disc of uniform field 
strength, so [3 for the warm gas is 100 times that of the cold 
(that is, the field is relatively weaker for the warm gas). An 
alternative approach to be explored in future calculations is 
to set up the cold gas in dense clumps such that [3 is uniform 
across the disc. 



3 RESULTS 

In this section we describe the structure of the disc for the 
single and multi-phase simulations. We ran these simulations 
for at least 300 Myr, by which time the structure of the spiral 
arms and morphology of the field were largely consistent 
with time. We discuss the morphology of the magnetic field, 
and the magnetic field strengths of the spiral arm and inter- 
arm regions in Sections 3.2-3.4. In Section 3.5 we compare 
the magnetic fields in our simulations with observations. 

3.1 Single phase simulations 

We performed a series of simulations of increasing field 
strength, taking a gas temperature of 100 K. In addition, 
we ran two calculations with warm (10^ K) gas. The param- 
eters for these models are listed in Table 1. Both the density 
and magnetic field strength increase during the simulations 
due to compression by shocks. The parameter /3 also changes 
and as described in Section 3.2.4, much of the gas in the disc 
exhibits lower values of /3 than indicated in Table 1. Further- 
more the calculations are scale free, in that the density can 
be increased, corresponding to a square root increase in the 
magnetic field strength, and the structure of the disc is the 
same. 

3.1.1 Spiral arm and inter-arm structure 

We show the column density for the calculations with 100 
K gas in Fig. 1. The panels (left to right, top to bottom) 
denote increasing magnetic field strength, or equivalent ly, 
decreasing j3 (as indicated in each panel). The degree of 
structure and the density of gas in the spiral arms clearly 
decrease as the magnetic field increases. For the model with 
the weakest field, which is essentially hydrodynamic, there is 
considerable inter-arm structure and the spiral arms them- 
selves show distinct clumps (Fig. 1, top left). The formation 
of spurs perpendicular to the arms and structure within the 
spiral arms has been previously desc ribed in terms of the or- 
bits of gas through the spiral shock (jPobbs k, Bonnellll2006l : 



iDobbs et al.l l2006l). Orbit crowding leads to gas of different 
velocities interacting in the shock and agglomerating into 
clumps. When (3 — 10, the substructure in the spiral arms is 
reduced, but spurs are still evident, and the density in the 
spiral arms is still relatively high. There are clumps along 
the spiral arms and substructure perpendicular to the spi- 
ral arms when /3 = 1 (magnetic pressure ^ gas pressure) 
but the density of these features is much less, and the inter- 
arm structure is much less distinct. For the case of a strong 
magnetic field, with /3 = 0.1 (Fig. 1, bottom right), there is 
a strong contrast to the models with a weak field (Fig. 1, 
top). There is little substructure and the spiral arms are 
much more continuous and smooth. The spiral arms are also 
much broader for the higher magnetic field strengths. 

Subsections of the disc for these results are shown Fig. 2, 
which more effectively highlight the inter-arm features and 
magnetic fields. The density contrast of the inter-arm struc- 
tures is clearer for the cases where the magnetic field is 
weaker than the gas pressure (top panels). These features 
are more distinct, and between the two inner arms, attain 
higher densities. For the (3 — ^.1 and f3 — 1 runs (stronger 
field), the structure of gas leaving the arms is much more 
continuous, hence the spiral arm appears broader. The re- 
duction in the strength of the shock also means that the 
inter- arm regions are more dense. Consequently the density 
of the inter-arm structure when [3 is 0.1 does not deviate 
significantly from the mean density in the disc and general 
dumpiness of the gas. By comparison there are relatively 
empty areas in the inter-arm regions in the essentially hy- 
drodynamic case and when = 10. 

The results from simulations using a gas temperature 
of 10^ K are shown in Fig 3. The top two panels show 
the case for [3 — 1. The shock is much weaker (in fact 
the gas barely shocks at all) compared to when the gas is 
cold, and the spiral arms are completely smooth and con- 
tinuous. The structure for the simulations where [3 — 100 
is not significantly diff erent from non-magnetic runs (see 
iDobbs Bonnelllbood l in which there is also no substruc- 
ture, but the shock is much weaker when [3—1 due to the 
magnetic field. 

3.1.2 Relation of the disc structure to the magnetic and 
thermal pressure 

The difference in structure between the magnetic and non- 
magnetic cases is largely attributable to the addition of mag- 
netic pressure (Eqn. 6). When the magnetic field is weak, 
this pressure is minimal and a strong shock occurs, similar 
to the hydrodynamic case. If B is large, the magnetic pres- 
sure reduces the strength of the shock and smoothes out the 
gas. This is similar in effect to increasing the thermal pres- 
sure, which is also found to reduce the structure in the disc 
(|Dobbs k Bonnelll [20Q6). 

For the case of a weak magnetic field, where va is 
low, the degree of structure depends primarily on the 
sound speed, as indica ted for the non-magnetic case in 
iDobbs Bonnelll d 2006h . When magnetic fields become im- 
portant, the strength of the shock depends on both Cs and 
VA- Although we show both warm and cold gas simulations 
with [3 — 1^ there is more structure in the case where the gas 
is cold since Cs and consequently va are lower. We performed 
a further lower resolution simulation (not shown) with cold 
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Figure 1. The spiral arm structure of the disc is shown for the single phase simulations with cold (100 K) gas after 250 Myr. The vectors 
indicate the magnetic field integrated through i.e. ^ Bxdz and f Bydz. The ratio of the thermal to magnetic pressure, f3, for each plot 
(also marked) is 10^ (top left), 10 (top right), 1 (bottom left) and 0.1 (bottom right). The degree of structure in the disc is reduced as 
the magnetic field strength increases. 



gas, taking P ^ 0.01, in order to give the same value of va 
as for the model with warm gas and p — 1. In this case 
magnetic fields dominate thermal pressure much more than 
expected from observations, but the structure is similar to 
that shown in Fig 3, top. Thus when magnetic fields are sig- 
nificant, increasing the magnetic pressure by a factor of 100 
has an equivalent effect to increasing the thermal pressure 
by the same degree. Assuming that the thermal and mag- 
netic pressures are equal (i.e. P ^ 1) we find the magnetic 
field has a strong effect in the simulations with cold gas, but 
does not remove all the substructure. 



3.1.3 Azimuthal profile and location of shock 

We plot the density against azimuth for the models with 
cold gas for P = 0.1, 1 and 10^, as well as the model with 
warm gas and /3 = 1, in Fig. 4. To determine the average 
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density we take gas in a ring situated at a radius of 7.5 kpc 
and width 200 pc which is divided into 200 sections over 
azimuth. The average density is calculated from the particles 
in each section. As expected, the density of the spiral arms 
increases as va decreases and the shock becomes stronger. 
We only show the average densities however - the maximum 
density is up to an order of magnitude higher than those 
displayed in Fig. 4. 

Also apparent in Fig. 4 is the shift in the location of the 
shock. The gas shocks further upstream for increasing va or 
magnetic field strength. Similarly the location o f the shock 
occurs further upstream at higher sound s peeds (|Slvz et alj 
l2003l : iGittins Clarkdliooi : lDobbsll2007l l. Thus increasing 
the magnetic field strength has a similar effect to increas- 
ing the thermal pressure, and the gas shocks earlier. For 
these simulations we find the gas generally shocks after the 
potential minimum. This has also been found in other re- 
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Figure 2. The spiral arm structure of the disc is shown for a 4 kpc x 4 kpc subsection of the disc shown in Figure 1, where the gas is 
cold. The ratio of the thermal to magnetic pressure, /3, for each plot is 10^ (top left), 10 (top right), 1 (bottom left), 0.1 (bottom right). 
The inter-arm structure becomes less distinct for higher magnetic field strengths, while for low field strengths, the magnetic field is much 
more disordered. 



cent results (iGomez Coxll2002l : IShettv Ostriked 120061 ) . 
the former using the same potential as in this paper. We 
note however that gas appears to s hock later for a 4 armed 
compared t o a 2 armed potential ([Gittins Clarkd I2OO4] : 
IPobbdlioOTl ) . The relative effect of including a magnetic field 
is to move the shock further upstream. For the case where 
[3 — 0.1, and the simulation with warm gas, we find that 
the shock is roughly coincident with the potential minimum 
where it would otherwise be downstrea m of the minimum 
without the presence of a magnetic field (jPobbsl 120071 ). 



3.2 Evolution of magnetic field for the single 
phase calculations 

Figures 1-3 also show the galactic magnetic fields in the 
single phase calculations. From Figs. 1 and 3, we see that 
on galactic scales, the field is strongly aligned with the spiral 



arms, and compressed in the spiral arms relative to the inter- 
arm regions. The exception is the case where the magnetic 
field is very weak (/3 = 10^), where the field is much more 
random. When there is a stronger field {(3 ^ 10), the field 
appears regular, and follows the large scale gas flow. For the 
model with warm gas and /3 = 1, (Fig. 3, top and middle), 
the very weak shock means there is only a small compression 
of the field in the spiral arms, and the field is less tightly 
constrained to the spiral arms. 

With a higher value of [3 — 100, and again a gas tem- 
perature of 10^ K, we find the morphology of the field is 
similar, but there is a stronger shock and a greater enhance- 
ment of the field in the spiral arms. In both simulations of 
warm gas, the field is very regular in the arm and inter-arm 
regions at both large and small scales. 

On smaller scales, (Fig. 2), we see that the magnetic 
field is not completely regular for the simulations with cold 
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2 4 

(radians) 

Figure 4. This figure shows the mean density plotted against 
azimuth for the simulations with cold gas, where ^ = 0.1 (blue), 1 
(green) and 10^ (red), and for the simulation with warm gas where 
f3 = 1 (magenta). The densities are calculated in the reference 
frame of the potential after 250 Myr. The azimuthal angle (j) is 
measured clockwise, in the direction of the fiow. The dotted lines 
indicate the minima of the potential. The density in the spiral 
arms decreases with increasing field strength, and the shock moves 
further upstream with respect to the potential minimum. 



arm regions, particularly when /3 = 1 and in some places 
the field is aligned in opposite directions to the large scale 
gas flow. As the field strength decreases (^^=10, Fig. 2, top 
right), the magnetic field becomes more disordered, even 
along the spiral arms. The field between the spiral arms is 
much weaker, and shows little coherence. For the essentially 
hydrodynamic case, {f3 = 10^) the field appears random ev- 
erywhere, although considerably enhanced along the spiral 
arms. At smaller size scales than shown in Fig. 2, the field 
is actually found to be more confined to the plane of the 
shock, where the gas is compressed, althoug h the field is 
aligned in both directions along the shock fc.f. lLaindll980l V 
If we zoom in to smaller size scales than indicated in Fig- 
ure 2, the field appears more disordered in the spiral arms 
for the simulations with cold gas. 



Figure 3. The spiral arm structure of the disc is shown for the 
single phase simulations with 10^ K gas after 250 Myr. The top 
panel shows the whole disc when (3 = 1 (i.e. magnetic and gas 
pressure equal), whilst the bottom 2 panels show a subsection 
when (3 = 1 (middle) and (3 = 100 (weaker field, bottom). The 
spiral arms are much smoother for the warm gas, but the spiral 
shock is very weak when the field is stronger (^ = 1). For both 
field strengths the field is strongly aligned with the gas fiow. 



gas, even with the strongest magnetic field {(3 — 0.1). We 
describe how the field disorder is generated in Section 3.3. 
In all cases except where (3 — 10^, the field is generally 
aligned along the spiral ams. However in the inter-arm re- 
gions, the field is more disordered. For [3 — 0.1 and 1 (Fig. 2, 
lower panels) , the field still largely follows the gas flow in the 
spiral arms, but there are significant deviations in the inter- 



3.2.1 Azimuthal profile 

Fig. 5 shows the average magnetic fleld strength plotted 
against azimuth after 200 Myr. The fleld strength is indi- 
cated for the simulations with cold gas where (3 — 0.1 and 
/3 = 10, and the simulation with warm gas and [3=1. We 
only consider the magnetic fleld in the plane of the disc 
here, i.e. \B\ = ^Bl+B^. We again determine the aver- 
age fleld from a ring of gas at 7.5 kpc as described in the 
previous section. For the simulations with 100 K gas, the 
fleld in the spiral arms is around 6 /iG when [3 = 0.1, and 
around 5 /iG when /3 = 10 (and 4 /xG when /3 = 1, not 
shown). The inter- arm magnetic fleld strengths are ^ 1/iG 
when [3^1. These are average values however - peak mag- 
netic fleld strengths are ^ 10/iG when /3 = 0.1 and [3 = 10, 
and ~ 20/iG when [3=1 for the cold gas. Although for the 
simulations with higher [3^ the magnetic fleld is weaker, the 
relative compression of the fleld in the spiral arms is greater. 
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Figure 5. The average magnetic field vs azimuth for the simula- 
tions with cold gas where /3 is 0.1 (red), 10 (blue) and with warm 
gas where (3=1 (green). The corresponding time is 200 Myr. 

For the simulation with warm gas, with a very weak shock, 
B increases by a factor of 2. This compares with factors of 5 
and 8 when /3 = 1 and 10 respectively, i.e. for progressively 
weaker fields. For the model with the very weak magnetic 
field, the mean magnetic field increases by up to 2 orders of 
magnitude in the spiral arms, although the mean strength 
in the spiral arms is only 0.1 /iG. The shift in the positions 
of the peaks in Fig. 5 again reflects the position of the shock 
with respect to the potential minimum. 

The relative amplitude of the magnetic field in the spiral 
arms is therefore related to the strength of the shock, as 
determined by Cs and va- We have already shown (Fig. 4, 
Section 3.1) that the density increase in the spiral arms is 
much greater when the field is weaker. Thus the increase in 
the magnetic field strength in the spiral arms will also be 
greatest when the magnetic fiel d is weak. W hen considering 
solutions to an oblique shock (|Priestl[l982l ). similar to the 
passage of gas through a spiral arm here, the ratio of the 
post and pre-shock magnetic field tends to the ratio of the 
post and pre-shock densities as va/cs 0. 

3.2.2 B vs z 

Recent observations from cosmic ray electrons suggest that 
magnetic fields exceeding 2 /xG are found up to at least 300 
pc above the mid-plane of our Galaxy, and are therefore 
stronger t han expect ed from the thermal pressure at these 
latitudes fc oxl liooHl ) . We find that the magnetic pressure 
largely follows the thermal pressure above and below the 
plane of the disc. Thus the magnetic field does not extend 
to such large scale heights, but we do not include supernovae 
or feedback in our simulations. 

3.2.3 B vs p 

We further checked for a relation between B and p, a 
property observed in the magnetic fields of galaxies and 
molecular clouds. For the gas in the spiral arms, we find 
B y9°"^^°'^ when P = 1 and the gas is cold (the exponent 



is slightly higher, ^ 0.75-0.8 for stronger shocks and lower, 
^ 0.65 for weaker shocks). This is comparable to the rela- 
tion expected for compression of a random magnetic field, in 
which case, by equating the magnetic and thermal pressure, 
B - y9^/^ For our simulations with cold gas, the magnetic 
field entering the shock has a significant random component 
assuming the gas has traversed at least one spiral shock. 
An alternative interpretation is that the relation is due to 
compression in the shock, with B scaling as p when B is 
perpendicular to the shock whilst B is independent of p 
when parallel to the shock. iGomez Cox, (2004 ) also find 
B oc /)°"^, which they relate to isotropic compression. Ob- 
servations suggest B oc p^'^^ for higher density gas, but for 
gas at lower densities, the exponent is much smaller (IValled 
Il995l ). 

3.2.4 The ratio of thermal to magnetic pressure, (3 

Although we take initial conditions of a uniform magnetic 
field and uniform density, at later times in our simulations, 
[3 takes a range of values about the initial value. However, 
much of the gas has a value of [3 lower than the average for 
the disc, particularly in the spiral arms. For example, where 
the initial /3 = 1, after 250 Myr, most of the gas in the 
spiral arms has /3 0.1, so the magnetic field dominates the 
gas pressure. In the inter-arm regions, the gas and magnetic 
pressure are more equal (so [3 ^ 1). Similarly for the case 
where the case where the initial /3 is 10, much of the gas in 
the spiral arms has (3 between 0.1 and 1. Since [3 oc p/B^ 
and we typically find B varies as p^"^, /3 tends to be lower 
for the dense gas in the spiral arms. Thus the initial values 
of (3 indicated in Table 1 (and likewise Table 2) only give a 
very approximate indication of [3 for the gas. 

3.3 The generation of a random magnetic field by 
spiral shocks 

Galactic magnetic fields are known to have both regular and 
random components (e.g. lBeck et al.|[l996h . Processes which 
generate turbulence in the ISM, i.e. stellar feedback, gravi- 
tational instabilities and MRI are expected to increase the 
random component of the field, but we focus on the effect 
of spiral shocks in this paper. 

We quantify the degree of order in the magnetic field 
by resolving the field into x and y components in the plane 
of the disc. The ordered and random components of B are 
calculated according to 

Bord = V< >2 + < By >2, 

Brand = V< [B^- < B^ + {By- < By >)2 > (13) 

SO that the total field is 

Btot = ^fBl^^+B^Zd- (14) 

However unlike the previous section, we volume average the 
magnetic field by integrating over the z direction of the disc. 
This is similar to Faraday rotation measurements and obser- 
vations from synchrotron emission. We take 200 points situ- 
ated evenly in azimuth round the disc at a radius of 7.5 kpc, 
and consider a cube of sides 100 pc centred on each point. 
The magnetic field in the Bx and By directions are found 
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Figure 6. The volume averaged random and ordered components 
of the magnetic field are plotted against azimuth for the models 
with cold gas for the strongest magnetic field ((3 = 0.1, top) and 
when the magnetic and thermal pressure are equal (^^ = 1, lower). 
The corresponding time is 200 Myr. The magnetic field is com- 
paratively more ordered when the field is stronger. 



by integrating over the z dimension at 100 x 100 points, and 
the averages taken over each set of 100^ points. 

Overall, we find that the magnetic field becomes more 
disordered as the strength of the shock increases, due to the 
increased velocity dispersion of gas in the shock. For a sim- 
ulation where there is no spiral potential, the field remains 
toroidal and of uniform strength for multiple orbits regard- 
less of Cs or VA- In our calculations with warm gas (where 
the shock is comparatively weak) the random component of 
the field is very small. For example, when /3 = 1 the random 
component is of order 1 per cent of the total field, whilst 
with (3 = 100 the random component is larger, but less than 
10% of the total field. However simulations with cold gas 
show a significant random component of the magnetic field, 
which increases as the magnetic field strength decreases, and 
dominates the ordered component for /3 ^ 10. 

Fig. 6 shows the random and ordered components of 
the field for simulations with cold gas, with /3 = 0.1 and 1 




200 
Time (Myr) 

Figure 7. The change in the ordered (solid) and random (dashed) 
components of the magnetic field, as a fraction of the total field 
are shown versus time for the models with cold gas for strong 
(f3 = 0.1) and moderate ([3 = 1) magnetic fields. The magnetic 
field becomes more disordered up to around 250 Myr as the gas 
passes through successive spiral shocks. 
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Figure 8. The random component of the magnetic field versus 
the velocity dispersion induced in gas in the shock after 250 Myr, 
for single phase cold (+) and warm (*) simulations. The random 
component generated is much larger when the velocity dispersion 
in the spiral arms is high. The points correspond to models A-F 
in Table 1. 



after 200 Myr. For the strongest field {(3 = 0.1, top), the 
ordered magnetic field exceeds the random field in both the 
spiral arms and inter- arm regions by a factor of 2-3. When 
/3 = 1 (producing a stronger shock) the random and or- 
dered components are similar in the inter-arm regions, whilst 
the ordered component is 1-2 times higher than the random 
component in the spiral arms. 

The degree of order in the magnetic field also changes 
with time, as shown in Fig. 7. The ordered component de- 
creases as the gas passes through multiple spiral shocks. The 
magnetic field becomes disordered in the centre of the disc 
first, as the evolution of this material is comparatively faster 
and the gas experiences a shock earlier. 
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Model Gas distribution Pcoid Pwarm B (/iG) 

G 50% cold, 50% warm gas 4 400 0.1 

H 50% cold, 50% warm gas 0.4 40 0.3 



Table 2. Table showing the parameters used for the two-phase 
calculations, where f3cold is the initial [3 for the cold gas, and 
likewise (3warm for the warm gas. 



3.4 Two-phase models 

The ISM is known to exhibit a multi-phase medium, so we 
now turn to models which include both warm and cold gas. 
Table 2 summarises the parameters for the two-phase mod- 
els. In this paper we do not include any heating or cooling be- 
tween the two phases. A resolution study (see Appendix A) 
shows that the morphology of interarm structure depends on 
resolution (spurs become more discontinuous and lose their 
integrity at higher resolution), but the typical spur spacing 
is independent of resolution. 



Note that the ratios of random to ordered fields depend 
somewhat on the averaging procedure. For example when 
using a mass weighted average for /3 = 1, the measured ran- 
dom component is a factor of 3-4 larger and thereby exceeds 
the ordered component. The random component of the field 
is reduced in the volume weighted average, since the field 
is most disordered in the high density regions at low scale 
heights. The mean field is also reduced by a factor of 1/2 
compared to Fig. 5 when using the volume weighted average. 

Recent simulations have shown that a velocity dis- 
persion can be induced in the ISM when c l umpy 
^as passes through a spiral shock (Bonnell et al.l l2006l : 
iDobbs k Bonnelll l2007al ). In IPobbs et al.l (|2006h , we de- 
scribe the dynamics of the shock as gas interacts and gains 
and loses angular momentum. For the MHD simulations pre- 
sented here, the shocks not only induce a velocity dispersion, 
but the velocities in the gas then lead to a disordered mag- 
netic field. For a low resistivity plasma (i.e. ideal MHD) 
the magnetic field is expected to follow the flow of the gas. 
Therefore the field will be become tangled where different 
velocities occur in a localised region of gas. We observe this 
in simulations where cold gas which has passed through the 
shock has a much more random field. We thus propose that 
spiral shocks generate and amplify a random component in 
the magnetic field. 

Fig. 8 shows the random component of the field plotted 
against the velocity dispersion induced in the shock. The 
dispersion was calculated for the velocity component in the 
plane of the disc, over subsections of a 200 pc width ring as 
was used to determine azimuthal profiles of the average den- 
sity and magnetic field. The velocity dispersions indicated 
in Fig. 8 are the average from the peak values in each spiral 
arm. For the simulations with warm gas, the gas entering the 
shock is fairly uniform, and a supersonic velocity dispersion 
is not induced in the shock (see Dobbs Bonne ll 2007a). By 
contrast, the velocity dispersion in the simulations with cold 
gas is highly supersonic, and consequently the random com- 
ponent of the magnetic field is much greater. As the mag- 
netic field increases for a given temperature (i.e. 100 K here), 
the magnetic pressure acts similarly to thermal pressure in 
both reducing the strength of the shock, and smoothing out 
inhomogeneities in the gas. Hence the velocity dispersion in 
the shock decreases and the field becomes more ordered. 

In some dense regions adjacent to the spiral arms 
(Fig. 2, lower panels), the field remains ordered, although a 
velocity dispersion is induced in the gas. This is presumably 
because there is low shear, so the gas and therefore magnetic 
field continue to follow the spiral arms. It is only on leaving 
these regions that the inhomogeneities in the magnetic field 
are sheared into more circular patterns 



3.4' i Spiral arm and inter-arm structure 

The results for the two phase calculations calculations are 
presented in Fig 9 with pcoid = 4 (left) and pcoid = 0.4 
(right). The top and bottom panels display the whole disc 
and a 4 X 4 kpc subsection respectively, the latter coincid- 
ing with the sections shown in Fig. 2. In both models the 
cold gas forms clumps, whilst the warm phase provides a 
background diffuse medium, similar to pre vious hydrody- 
namic simulations (iDobbs &: Bonnell l2007bl V The densities 
of the spiral arms and inter-arm clumps are reduced for the 
stronger fields (lower /?), whilst for weaker fields (higher 
P), the clumps are more coherent and form longer, more 
spur-like structures. From the left hand panels in Fig. 9, 
we see that the inter-arm gas is similar to the 'rungs of a 
ladder' between the spiral arms, with dense clumps situated 
along them. This structure is persistent with time, as demon- 
strated in Fig. 10, which shows the disc in the f3coid = 4 run 
after 600 Myr. The structure in this case is much more sim- 
ilar to M51, the archetypal galaxy for comparing inter-arm 
features, which has long thin spurs dotted with star forming 
regions. 

Whilst the magnetic pressure tends to weaken the spi- 
ral shock and smooth out structure, the presence of a warm 
phase counteracts this to some extent. For the case where 
Pcoid = 4, the inter-arm clumps are more dense than the 
corresponding time in the single phase model where /3 = 10 
(Fig. 2, top right). We have previously noted in hydro- 
dynamical simulations that the warm component provides 
a pressure which confine s the cold gas to higher densities 
(jDobbs k Bonnellll2007bl ). Overall we find from the simula- 
tions in this paper, and lower resolution simulations, that 
if Pcoid ^ 0.1, the cold and warm phases do not separate 
(since even in single phase simulations the magnetic pres- 
sure smoothes out any structure in the cold gas), but tend 
to remain much more randomly distributed. With f3coid ^ 1, 
more coherent features develop in the cold gas, but for in- 
termediate values only smaller clumps tend to form. 

Thus in the results presented in this paper, the cold and 
warm gas are essentially separate, i.e. the cold gas is situated 
in clumps embedded in the warm phase. This is because the 
increase in density of the cold gas due to the spiral shock is 
much larger than that of the warm phase. Orbit crowding 
in the spiral shock then leads to the agglomeration of these 
clumps into larger structures, though this process is reduced 
by magnetic pressure, and is more evident when self gravity 
is included (Dobbs et. al., in preperation) . This occurs re- 
gardless of whether the simulations are single or two phase. 
The densities of the cold gas are higher when a warm phase 
is present (particularly in the inter-arm regions), but the 
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Figure 9. The column density is shown for the two-phase simulations after 250 Myr, for the whole disc (top) and a 4 x 4 kpc subsection 
(bottom). The left hand panels show the case where Pcold = 4 and the right hand panels where f3coid = 0-4. Both the cold and warm 
phases are shown in the plots, but we show them separately for the case where f^coid = 4 in Fig. 12. There is more structure in the cold 
gas when the magnetic field is weaker {^cold = 4). The vectors show the magnetic field smoothed over a particular grid size. There is 
more detailed structure on smaller scales, particularly in the spiral arms which are better resolved. 



structure of the cold gas is similar (e.g. Fig. 9 top and Fig. 2 
bottom left). 



3.4-2 Evolution of magnetic field 

Similarly to the single phase results, the magnetic field ap- 
pears regular on galactic scales, but more disordered in the 
inter-arm regions on smaller scales (Fig. 9). There is again 
a more random magnetic field in the model with a higher (3 
(Fig. 9, left, compared to Fig 9, right). In Fig. 11 we plot 
the cold (top) and warm (lower) components separately for 
the model with f3coid = 4. As expected from the single phase 
results, the magnetic field of the warm gas is more regular 
than that of the cold, especially in the inter-arm regions. 
However the magnetic field of the warm gas is clearly more 



random than in entirely single phase simulations of warm 
gas (compare to Fig. 3). 

We show the volume averaged ordered and disordered 
components of the magnetic field for the two-phase simula- 
tion with Pcoid = 4 explicitly in Fig. 12. For the warm gas, 
the disordered and ordered components of the field are of 
comparable strengths. This is in contrast to the single phase 
simulations of warm gas, where the disordered component 
was much smaller than the total field. Thus the magnetic 
field in the warm gas is modified by the distortions of the 
field in the cold gas and the velocity dispersion induced in 
the cold phase. 

A caveat to these results is that since we took initial 
conditions consisting of a uniform density disc of uniform 
magnetic field strength, (3 for the warm gas is 100 times 
higher than that of the cold (an alternative would be to 
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Figure 10. The column densities and magnetic field for the twc 
phase simulation with f3coid = 4 are shown after 600 Myr. 

start with an initially clumpy distribution). As mentione 
in Section 2.1, p is expected to be similar for the cold an 
warm components of the ISM. However, a smaller value ( 
(3 for the warm gas is unlikely to change the structure of tli 
disc significantly in the two-phase results since this primaril 
depends on the cold gas, and the warm phase merely aci 
as an extra pressure. Although the magnetic field may t 
more ordered with a lower /3, we find the field is much moi 
disordered when p for the warm gas is 40 in the two-phae 
results compared to /3 = 100 in the single phase result: 
Thus it seems likely that a highly irregular field in the col 
gas will have some effect on the field in the warm gas eve 
if (3 is lower for the warm gas. 

3.5 Comparison of magnetic field with 
observations 

3.5.1 Warn Neutral/ Ionized Medium 

The principal technique for measuring magnetic fields in the 
warm ISM is from radio measurements of diffuse synchrotron 
emission. Synchrotron emission arises from the interaction 
of cosmic rays with warm gas and has enabled mapping of 
magn etic fields in several external galaxies such as M51 and 
M81 (|Beckll2007l ). An estimate of the distribution of cos- 
mic ray electrons is also required in order to determine the 
magnetic field strength. For A = 6 cm, Faraday rotation 
of the emission is negligible, and radio polarisation provides 
the strength and orientation of the regular component of the 
field. In conjunction with Faraday rotation measurements at 
several wavelengths, the direction of the field can be mea- 
sured. 

Observations from synchrotron emission are generally 
associated with the warm component of the ISM. The to- 
tal (random + ordered) magnetic field is typically ^ 6/iG, 
although magnetic field strengths are higher in the spi- 
ral ar ms of grand desig n galaxies, with strengths of lO's 
of //G (|Beck et al.lll996h . Our simulated field strengths for 
the warm gas tend to be lower by a factor of 3 to 4, al- 




-6 



Figure 11. The column densities and magnetic field for the cold 
(top) and warm (gas) are shown separately where (3coid = 4. This 
4x4 kpc subsection of the disc is the same as the bottom right 
panel of Fig. 9. The magnetic field is much more disordered in 
the cold gas compared to the warm component. 

though the maximum field strengths are much more compa- 
rable to the synchrotron measurements. The mean strengths 
are closer to Faraday rotation measurements of pulsars in 
the Galaxy, which indicate lower field strengths (e. g. 3/i G, 
(IHan, Manches ter, Lvne, Qiao van Stratenlliooeh ). How- 
ever the magnetic field scales with density. In providing val- 
ues for the magnetic field, we took a surface density com- 
parable to the solar radius, wh ere the magnetic field is rel- 
atively weak (|Han et al.l 120061 ) . The random and ordered 
components of the field are thought to be of similar magni- 
tude. This is the case for the two phase simulations, whilst 
the random component is minimal in simulations with just 
warm gas. 
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Figure 13. A synthetic polarisation map for the two phase simu- 
lation where [Scold = 4. The arrows (B— vectors) are constructed 
from the Stokes Q and U parameters and the contours show the 
total synchrotron intensity (Stokes I parameter) smoothed to a 
beam size of 0.7 kpc (ie. /ibeam = 0-35 kpc). 



into the inter-arm regions. By contrast in our simulations, 
the enhancement of the magnetic field is strongly confined 
to the regions corresponding to the spiral shock. Some 
observations of galaxies even find that the magnetic spiral 
arms do not coincide with the optical spiral arms where 
star formation is occurring (Beck Hoernes 1996). This 
phenomenon is difficult to explain purely by the dynamics 
of spiral shocks, as is evident from our models. 



Figure 12. The volume averaged ordered and disordered com- 
ponents of the magnetic field are plotted for the warm and cold 
gas in the two-phase simulation with Pcoid — 4 for the cold gas. 
The top panel shows the magnetic field for the cold gas, whilst 
the lower panel shows the magnetic field for the warm gas. 

3.5.2 Synthetic observations 

In order to make a meaningful comparison of our results 
with observations of magnetic fields in spiral galaxies, we 
have computed a synchrotron polarisation map of one of 
our simulations. A detailed description of how the polarised 
emission is calculated is described in Appendix B. 

The synthetic polarisation map is shown in Fig. 13 
for the two phase simulation where pcoid = 4. The po- 
larised emission is enhanced in the spiral arms, whilst the 
magnetic field is generally aligned with the spiral arms 
comparable to radio- intensity maps of M 51 (Neininger 'l99^ 
Patrikeev, Fletcher, Stepanov, Beck, B erkhuiisen. Frick & £ 
20061 : lBeckll20Q7l ). In our simulations, the spiral shock cor- 
responds to the dust lanes of actual galaxies, generally 
coincident with trailing edge of the optical spiral arms. 
However the most notable difference from our simulations 
is that the polarized intensity i n the observations of M51, 
and similarly M83 (|Beckl I2OO7I ) forms much wider spiral 
arms, with contours of the emissivity spread over a few kpc 



3.5.3 Cold Neutral Medium 

Unfortunately there are no corresponding observations of 
the magnetic field in cold HI on galactic scales. Zeeman split- 
ting has recently been used to determine magnetic fields in 
cold HI structures in the Arecibo survey (Heiles Trolandl 
I2OO5I ). but this method requires strong magnetic fields. Fur- 
thermore the observations do not provide spatial informa- 
tion on the magnetic field, rather a median field strength 
of 6/i G. The cold gas in the spiral arms of our simulations 
has a field strength of around 2/i G in the two phase sim- 
ulations and with (3—1. The relatively low magnetic field 
strengths in our simulations suggest that it would be worth 
setting up simulations where the cold gas is already situated 
in dense clumps, and establishing whether this structure is 
maintained for higher field strengths. We find a compara- 
tively greater increase in the magnetic field strength in the 
shocks for the cold gas than warm gas, typically by a factor 

^ 5 -10. 

I Li et al.l (|2006l ) calculate order parameters (i.e. the ra- 
tio of the ordered to total field) of ~ 0.2 — 0.4 in several 
~ 100 pc size GMCs, where the configuration of the mag- 
netic fields in the GMCs is thought to be a consequence of 
the external magnetic field rather than internal sources. We 
expect more detailed calculations to show whether GMCs in 
computational models exhibit a similar degree of order. 
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4 CONCLUSION 

We have performed simulations of galactic discs subject to 
a spiral potential with a range of field strengths. The main 
results we have discussed are 1) the reduction in structure 
across the disc as the magnetic field strength increases, and 
2) the possibility of spiral shocks inducing an irregular mag- 
netic field in the ISM. 

As the strength of the magnetic field increases, the 
strength of the spiral shocks and therefore density of the 
spiral a rms are reduced. This is as expected from the anal- 
ysis of Roberts Yuaiil (|l970l V Consequently the forma- 
tion of spurs is increasingly suppressed for higher magnetic 
field strength. The spiral arms themselves are more con- 
tinuous and clumps along the spiral arms are less dense. 
This supports prev ious 2D results (She ttv Ostr iker 2006; 
iTanaka et al.|[2Q0^ ) which find that both Kelvin-Helmholtz 
and gravitational instabilities perpendicular to the magnetic 
field are reduced by stronger magnetic fields (although we 
do not relate the structure in our simulations to these in- 
stabilities, our overall conclusions agree). We find that the 
addition of the magnetic field is similar in effect to an in- 
crease in thermal pressure, in that both provide a pressure 
which oppose the formation of structure, and smooth out 
the gas. 

Nonetheless, we still find significant inter-arm struc- 
ture with the presence of a magnetic field , unlike the results 
(those which are non-self gravitating) of IShettv Ostriker! 
(|2006l ). This structure is present in the cold component of 
the ISM, which has generally not been included in previous 
simulations. Inhomogeneities present in the initial random 
distribution of gas become amplified by spiral shocks. With 
warm gas or strong magnetic fields, these inhomogeneities 
are smoothed out by the pressure. For our calculations with 
cold gas, we find that magnetic fields only prevent substruc- 
ture when the ratio of gas to magnetic pressure {(3) is < 0.1. 
For /3 > 1, spurs perpendicular to the arm still form in the 
cold gas. The additional pressure provided by warm gas in 
the two-phase results increases the longevity of structure in 
the inter-arm regions, even when the magnetic field domi- 
nates for the cold gas. 

The main difference between ours and previous results 
is that we use cold gas, although we also use a weaker mag- 
netic field than lShettv Ostriker! (|2006h . We find structure 
is present for values of f3 in the cold gas similar to observa- 
tions, but the volume averaged magnetic field strengths in 
our two phase models are typically < 1/iG, lower than ob- 
servations. This discrepancy arises because we start with an 
initially uniform distribution of cold gas. To obtain a similar 
[3 for both the cold and warm gas in our two phase models, 
we would need the cold gas to be 100 times denser than the 
warm gas, i.e. this just indicates that structure is already 
present in the cold gas. If we u sed higher field strength s, 
we would find, like the results of Shett v Ostri ker (2006), 
very little structure emerges. In this case, (3 for the warm gas 
would be closer to observations, but the magnetic pressure 
would dominate the thermal pressure by a factor of 100 or 
so in the cold gas, inconsistent with observations. Thus we 
either need to distribute the cold gas initially in clumps, or 
ideally include a much more detailed thermal treatment of 
the ISM, in order to obtain magnetic field strengths and f3 
consistent with observations. 



We find that whilst the ISM appears highly structured 
in observations, we would expect a higher degree of struc- 
ture with relatively weak magnetic fields. For instance, the 
two phase model where pcoid = 4 retains much more struc- 
ture typical of grand design galaxies compared to the case 
where f3coid = 0.4. Current observations s uggest that mag- 
netic pressure exceeds thermal pressure (iHeiles k Trolandl 
I2OO5I I . We however note that /3 exhibits a range of values in 
our simulations, and (3 tends to be take smaller values in the 
spiral arms and dense gas (for a given temperature), which 
are more likely to correspond with observations. Again a 
more complete treatment of the ISM in future work, and 
further observations of the CNM will allow a better com- 
parison. 

In our simulations, the magnetic field is compressed by 
the spiral shocks, a s expected from a straightforward a naly- 
sis of MHD shocks ([Roberts Yuanlll970l : IPriestll 19821 ). The 
relative increase in the magnetic field strength is greater 
where the shock is stronger. However the most intriguing 
result from our simulations is the possibility that spiral 
shocks generate an irregular magnetic field. This process 
has not been identified in previous simulations, which we 
attribute to the fact that they have not included a cold 
phase. Galaxies are known to contain a random compo- 
nent of the field, but it is usually supposed that this is 
due to supernovae and/or feedback from stars. We there- 
fore postulate that spiral shocks are important in gen- 
erating disorder in the magnetic field, whilst simultane- 
ously inducing a velocity dispersion and density struc- 
ture in the gas (|BonnelL Dobbs, Robitaille Pringlell2006l : 



iKim, Kim Ostriker! I2OO6I : iDobbs et al.ll2006l ). The degree 
of order in the disc, similar to the presence of arm/inter-arm 
structure is related to the strength of the shock. When the 
gas is cold and the magnetic field is weak, the shock is much 
stronger, and a higher velocity dispersion is induced in the 
gas. Consequently the magnetic field lines follow the gas and 
the field becomes tangled. In simulations with warm gas, we 
find the magnetic field is almost entirely regular, although 
with a stronger shock, the field may become more irregular. 
However in our two-phase results, the velocities and irregu- 
larities in the magnetic field of the cold gas induces compar- 
atively more disorder in the warm gas than the single phase 
calculations. Although there are no measurements of regular 
versus disordered components of the magnetic field in cold 
HI, we would predict from our results that the field in the 
cold gas will be more disordered than that of warm HI. 
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APPENDIX A: RESOLUTION TESTS AND 
SPACING OF SPURS 

In order to demonstrate that the length scales of structure 
in our simulations is not dependent on resolution, we per- 
form a study of the two phase simulations with pcoid = 4 at 
different resolutions. Fig. Al, shows simulations with 1 (left) 
and 8 (right) million particles. With higher resolution, the 
structure tends to become more disjointed, and therefore it 
is harder to recognise large scale features. However the main 
spurs which extend across the inter-arm regions tend to be 
similar at different resolutions. In fact the structure of the 
cold gas in the two phase simulations is essentially the same 
as an equivalent simulation with only cold gas. 

The large scale structure and spacing between spurs is 
determined by the dynamics of the spiral shock, and there- 
fore should not be affected by the numerical resolution. In- 
stead we find that the spacing between spurs is dependent 
on the strength of the spiral shock, and is thus related to 
the epicyclic radius, and the thermal and magnetic pressure. 
Given that the primary purpose of the present paper is to 
examine the effect of including a magnetic field, we intend to 
discuss this in an upcoming paper (Dobbs, in prep.) which 
investigates a range of parameters, including the strength 
of the potential, the addition of self gravity, and in these 
calculations, the spacing between spurs is clearer. 

We have also performed calculations where the cold gas 
is initially distributed in cold clumps, rather than randomly. 
In this case, the initial structure is somewhat different be- 
cause the clumps are sheared into thin dense filaments. How- 
ever as the gas passes through the spiral shocks, this struc- 
ture tends to be disrupted, and the spacing between spurs 
is similar to that apparent in Fig. Al. 



APPENDIX B: CALCULATION OF 
POLARISED SYNCHROTRON EMISSION 

The method for computing the sample po- 
larised synchrotron emission map shown in 
Section 3.5.2 is describ ed here, follow ing 
ISokoloff, Bykov, Shukurov . Berkh uiisen, Beck Poezdl 
(|l998l ) and lUrbanik. Elstner k Beckl (|l99/t ). 

We adopt a reference frame such that the galaxy is ex- 
actly face-on to the observer, such that z is the coordinate 
along the line of sight and x,y are the coordinates both 
in the plane of the sky and also in the galactic plane. The 
direction (B vectors) and intensity of synchrotron polarisa- 
tion are determined observationally by the Stokes Q, U and 



I parameters (respectively), given by 
Q = Pi J i(;(r)s(r) cosxdV, 
Pi J w{r)s{r) sin x^V, 
I ^ p, j w{v)e{v)dY 



U 



(Bl) 
(B2) 
(B3) 



where the integration is over the cylindrical column consist- 
ing of an assumed beam profile w(x,y) along the line of sight 
(z), e is the local synchrotron emissivity, pi is the intrinsic 
degree of polarisation and the angle x is related to the local 
magnetic field according to 



tan~\By/Ba:) + KX^ J B^n(r,z)dz. 



(B4) 



The second term in the above expression introduces the ef- 
fect of Faraday rotation: a wavelength- dependent (oc A^) 
change of the polarisation angle dependent on the strength 
of the magnetic field parallel to the line of sight (Bz) and 
the number density of thermal electrons, n(r, z) (the param- 
eter K is a constant given hy K = 0.81 rad m~^ cm^/iG~^ 
pc~^). Since Faraday rotation is wavelength- dependent it a) 
becomes negligible at short wavelengths (ie. for observations 
made shortw ard of A 6cm for typical ISM magnetic fields 
(jBeckl l2Q07l ) and b) can be corrected for by observations 
made at multiple wavelengths. Thus we show our results 
assuming no Faraday rotation. 

The synchrotron emissivity is given by: 



.{r,z) 



(B5) 



where a is the synchrotron spectral index (Ii, oc z^~"), 
B± = B^ + B^ is the magnetic field in the plane of the 
sky and ncr{r,z) is the local number density of cosmic ray 
electrons. The intrinsic degree of polarisation, pi is related 
to the spectral index hy pi = (a + l)/(a + 5/3), which for 
our choice of a = 0.8 gives pi = 0.73. 

The electron cosmic ray distribution is assumed to be 
of the form 



: exp[-r/rc 



z/ Zcrit] 



(B6) 



where the distribution is assumed to fall off with galactic 
radius assuming rcrit = 13.5 kpc and with galactic scale 
height using Zcrit = 2.5 kpc (|G6mez Coxl [20041 : iFerrierd 
[1998). 

In order to construct the polarisation map efficiently 
from our SPH simulations we use a beam profile w b ased on 
the usual SPH cubic spline kernel (|Monaghanlll992l '). 



Wiq) 




0^q<l; 
1 ^ Q < 2; 
q>2. 



(B7) 



where q — ^ (x^ + + z'^)/h. This kernel is approximately 
Gaussian in shape but with compact support of radius 2/i, 
where h is the smoothing length). We perform the integra- 
tion by interpolating the quantity to be integrated at a given 
position from nearby SPH particles, ie. 



{A{x,y,z)) = V^^^-mj(x,2/,2;,/i, 



(B8) 



where m^, pj and hj are the mass, mass density and smooth- 
ing length of particle j, each of which are known from the 
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Figure Al. The column density is shown at different resoutions for the two-phase simulations with [Scold = 4 after 250 Myr. The 
reolution is 1 million particles (left) and 8 million particles (right) for the whole disc (top) and a 4 x 4 kpc subsection (bottom). At 
higher resolution, more detailed structure emerges, and the spurs become more fragmented. However the large scale features are similar 
at the two different resolutions. 



simulation. For example in the case of the Stokes Q param- 
eter, we have A = e{x,y,z) cosx (similarly for the U and / 
calculations) . 

The use of the summation interpolant (equation IB8|) 
takes care of the x,y dependence of the integral (having re- 
placed it with a summation over contributing particles) . The 
line-of-sight (z) integration in ()B1|) - (|B3|) is then performed 
by integrating ()B8|) analytically, giving 

J (A)dz = Y,^A, J W^j(x,y,z,h)A^, (B9) 

such that the z integration simply requires the integral of 
the kernel through one spatial dimension. This can be pre- 
computed efficiently in a table, giving, in effect 

j{A)dz = Y,j^A,Y,Ax,y,h), (BIO) 

where Y \s, W integrated along one spatial dimension. To 
approximate a beam profile w(x,y) = Y(x,y) with a fixed 
resolution length (rather than the varying resolution length 
used in the SPH simulation), we set the smoothing length 
used in the interpolation to /i = maj^{hj , hteam) - In prac- 



tice our resolution is such that hbeam » hj (that is the 
simulation is much better resolved than the observations). 

The net result of the procedure given above is that the 
construction of the polarisation map is reduced to a loop 
over all the particles in the simulation, where for each par- 
ticle a contribution is added to all pixels in the image lying 
within 2h. 

The additional inclusion of Faraday rotation would add 
a slight complication to the calculation in that the angle 
X also depends on a line-of-sight integration. This would 
require sorting the particles in z prior to computing the sum 
()BlQp , then computing the sum from back to front with x 
stored and incrementally updated for each pixel. 

The Fortran routine for performing this interpolation 
has been incorporated into SPLASH and is available on re- 
quest. 
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